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ABSTRACT 

All present stiffened cylindrical shell design formulas 
for the case of external hydrostatic pressure were surveyed. 
This also included present design practices for allowances 
due to pressure hull imperfections, and actual test data 
when available. Formulations for all. three basic hull 
failure modes were then selected, first for accuracy, and 
secondly for compatibility with both elastic-perf ectly 
plastic and strain-hardening metals, whenever possible. 

The formulas were then inserted in a logical flow 
pattern to design elastic-perf ectly plastic scantlings .for 
failure at the most efficient collapse mode. The process 
was programmed in FORTRAN IV and designed to iterate, vary in 
several scantling parameters systematically. The 
"optimum" design, based on a simple hull weight/buoyancy 



ratio was selected at the completion of the run. 
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Program inputs are: collapse depth, hull diameter, 

hull length, internal bulkhead spacing (specified or 
unspecified), framing (internal or external), and metal 
properties. Outputs are: shell thickness, typical frame 

spacing, typical frame size, heavy frame (bulkhead) spacing, 
and heavy frame size for each design, and an optimum design 
designation. Simple directions are given for conversion 
of the program to one compatible with strain-hardening metal 

Thesis Supervisor: J. Harvey Evans 



Title; Professor of Naval Architecture 
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LISTS CF SYIi^OLS 
USED COMiOMLY IN THE TEXT 



( i--v/2') a 

shorthand for express ion :used in ELI.Gh 

displacement of pre-specif ied length of hull, tons 



hull mean diameter (to shell mid-fiber), in. 
Young’s Modulus, psi 
secant modulus, psi 
tangent modulus, psi 

unsupported length of plating between frames, in. 
applied hydrostatic pressure, psi 
hull mean radius, in. 
shell thickness, in. 

weight of pre-specif icd length of hull, tons 
measure of beam-column effect (see chapter 2) 
shell fexibility parameter (see chapter 10) 

4 I r— 7- 

"thinness" factor, /( L/Pn ) 2 / Oy 

v( t/PnT3 \/ E* 



Poisson's ratio 

circumferential bending stress, midbay, psi 



f circumferential] 
( longitudianal •/ 

yield stress, psi 



membrane 



stress , 



nidbay , 



psi 



USED IN PROGRAM 



sane as (1) above. 
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20. AE: 

21. AF: 

22. AFH: 

23. ALFA: 

24. B: 

25. BETA: 

26. BH: 

2?. BHETA: 

28. B3 : 

29. CC: 

30. CCOP : 

31 . D1SP : 

32 . BH : 

33. DK: 

34. BLH: 



effective cross-sectional area of typical 
frame, sq. in. 

actual cross-sectional area of typical frame, 
sq. in. 

actual cross-sectional area of heavy frame, 
sq . in. 

AE/(FB*T), used in definition of A. 
typical frame web thickness, in. 

B/F3, used in definition of A. 
heavy frame web thickness, in. 

Von Sanden-Gunther variable vised in RFIDS 
(included so as not to be confused with BETA) 
Bulkhead spacing (i.e., heavy frame spacing), in, 
typical frame dimension parameter (see chapter 10 
and Appendix A) 

typical frame moment of inertia parameter derived 
from CC. 

same as (2) above 
same as (3) above 

distance from shell mid, -fiber to combined 
centroid of typical frame and effective length of 
shell plating, in. 
same as. (33). for heavy frames. 

C, ETAM : same as (4), (5) and (6) above. 
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40. 

41. 

42. 

43. 

44. 

45. 
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47. 

48. 

49. 

50. 
51 • 

52. 

53. 



EX: effective moment of inertia, of typical frame and 

effective length of shell plding , in. 

EL: effective length of shell plating, in., for typical 

frames. 

ELh: effective length of shell plating, in,, for heavy 
frames. 

Fl , F2, F3, F4: Salerno-Pulos "F-functions" used in 

i 

computing shell stresses. 

FC: distance from shell mid-fiber to centroid of 

typical frame, in. 

FCH: same as (40), for heavy frames. 

FD: typical frame depth, in. 

4 

Fi : typical frame moment of inertia, in. 

FIH: typical heavy frame moment of Inertia, in. 

F8: typical frame spacing, in. 

FW: typical flange vji-dth, in. 

GAI-iA. : same as (12) above, 

Gi.U: Poisson’s ratio 
KULNTH: hull length, in. 

UP : number of typical frames 

PC: design collapse pressure, psi 

?CG : general instability collapse pressure, psi, 
reduced for imperfections. 

PCGE1 : general instability elastic collapse pressure, 

psi 
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55 . 
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58. 

59. 

60. 
61 . 

62. 
6 3. 

64. 

65. 

66 . 
•67. 
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70. 

71. 

72. 
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74. 



PCLSl : axisymmetric yield collapse pressure, psi 

PC R: asymmetric buckling pressure, psi, reduced for 

imperfect ions. 

PCF1S or PCH21 : asymmetric elastic buckling pressure, 

psi 

PEL: axisymmetric elastic buckling pressure, psi 

PO: operating depth pressure, psi 

PRE: P. e , chapter 3 

PRI : P^ , chapter 3 

RC: hull radius to centroid of heavy frane-ef f ective 

length of shell plate combination, in. 

RCG: same as (6l), with typical frames 

RF: hull radius to centroid of typical frame, in. 

RHO: material density, lb/ in . 3 

KM: hull mean radius (to raid-fiber of shell), in. 

RO: hull outer radius (to outer fiber of shell), in. 

S1GY : yield stress, psi 
SL: same as (7) above 

T: same as (10) above 

TF: typical frame flange thickness, in. 

THETA: same as (13) « above 

VF: volume of typical ring frames; in^ 

VFK: volume of heavy ring frames, in^ 

V/D: weight/displacenent ratio 
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optimum WD 
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WT : 


same as (11) above 
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input parameter specifying frame location 
(i.c., 1 « 0 rr) internal , anything else zz) external ) 
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INTRODUCTION 

The amount of literature concerning the collapse of ring- 
stiffened cylindrical shells under hydrostatic pressure 
accumulated in the last fifty years is voluminous. This is 
understandable; the subject is very involved. To this day, 
an exact solution for all aspects of the problem does not 
exist. Goo d solutions for the different failure modes do 
exist, though, modified in varying amounts by empirical data. 
No attempt will be made here to list or summarize this 
knowledge. Many have already done this, and the finest 
review to date in this author's knowledge has been done by 
J. G. Pulos 1 ^ for the Navy's former David Taylor Model Basin. 

It appeared that one should be able to integrate this 
albeit incomplete, yet extensive knowledge with the use of 
present generation computer science. Hand calculations for 
only one geometrical combination of shell thickness, frame 
size, frame spacing, hull diameter, hull length, etc. are 
notoriously laborious, even for only one mode of failure. 
Submarine design processes using the hand technique would 
achieve adequate structures, but v?ith little or no idea if 
anything better existed. ' Optimization, vrith the exception of 
a few combinations tried at a great cost in time, was out of 
the question: while similar submarines could be designed on 

past knowledge, different hull geometries or deeper operating 
depths meant a great deal of time and work. It was at the 
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suggestion of Professor Evans that the development of a 
computerized design optimization v?as undertaken. The basic 
design equations were there. The computing tools were 
readily available. All that needed to be done was an 
integration of the two into a practical, useable, and most of 
all, reliable (in terms of latest empirical data, if 
necessary) program. 

Shortly after the start of the project, it was 

p p 

discovered that a smiliar program had just been completed^ 
for the Naval Ships Research and Development Center (NSRDC). 
It is hoped that by using some different approaches and 
techniques, this program could be a valuable tool to use in 
conjunction with reference 22. 

The description of the program development will be done 
by subprograms, each building on the other, and ending with 
the optimization scheme of the main program. While the 
program developed can be used only with elastic-perf ectly 
plastic, isotropic materials (e.g., HY-30 steel), it can 
easily be modified so as to be applicable also to strain- 
hardening metals, such as HY-150 steel. Further discussion 
on this will follow’ later. Included also in the thesis will 
be various data obtained using the program or portions 
thereof in parametric-type studies. 
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It is recommended that while reading through the various 
chapters on subprograms, reference be made to figure 23 
(chapter 10), which is a simplified main flow diagram for 
the entire optimization. 
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CHAPTER Oi'iE 

BACKGROUND INFORMATION AND TERMINOLOGY 
The entire hull design program is based on the three 
fundamental failure nodes for ring-stif f ened cylindrical 
shells. Therefore, in order that terminology remain 
consistently clear throughout the discussion, a brief 
description and catagcrization of these modes follows. Refer 
to figure one for pictorial representations. 

A . AXI SYMMETRIC FAILURES BETWFBi: STJ.FFE KEPS 

All axisyrametric failures, whether elastic, plastic, 
or some mixture thereof, are characterized by one or more 
accordian-like pleats, or circumferential ripples between 
ring frames. For true axi symmetric failures, the stiffeners 
remain undeformed. 

1, AXISYMKETRIC YIELD FAILURE. This type of failure 
occurs only with elastic-perf ectly plastic 
(plateau-type stress-strain curve) materials. It 
is regarded, then, as almost totally a yield-type 
failure, although it is initiated partially by 
instability phenomena. 

2. AXI SYKkET R I C •113 1 . AST I C FAIL URE. Inelastic 
failure is also a failure above the purely 
elastic range, but in strain-hardening materials. 
Thus, the failure is in the range where Young’s 
modulus- varies , and can intuitively be 
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considered a kind of combination elastic-plastic 
failure. Often this failure can be at a lower- 
pressure than a pure buckling (elastic) failure, 
due to certain combinations of lov: modulus and 
hull geometry. The pure yield failure (A.l.) 
cannot occur, .in strain- hardening materials. 

3 * AXI S YMMETRIC ELASTIC FAILURE . This is 

axi symmetric failure in the linearly elastic 
range. Theoretically, this could occur, given 
proper geome.try, in either of the tv?o above types 
of materials. Generally speaking, the required 
geometry is one of a thin shell relative to the 
hull diameter and depth. In reality, this 
failure is, at this writing, a mathematical 
phenomenon only. Other hull failure modes occur 
first, so this mode has never been achieved in 
actual testing. It is valuable, however, in 
determining effects of geometrical defects on 



collapse pressure. 



B. ASYMMETRIC. .(LOBAR) BUCKLING FAILURES B ETWEEN 
STIFFENERS : 

All lobar buckling failures are characterized by- 
lobes of buckling distributed partly or completely around the 
shell circumference. The frames remain intact. 



1. ELASTIC OR- INELASTIC MODES : 



is primarily 



buckling failure 



The collapse pressure is 
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dictated by the modulus (elastic or inelastic) 
and hull geometry, 

C . GEMER tl, INSTABILITY ( SHELL AID F ? AKS ) FAILURE . 

General instability failures are characterized by 
failure of both frames and shell simultaneously, sometimes 
extending the entire length of the cylinder. 

1. ELASTIC Oh I N BLAST I C MODES . These are also 
, buckling failures, except that lobes extend 
longitudinally as well as circumferentially. 

It is normally assumed that only a one-half wave 
extends between heavy frames (bulkheads), since 
these heavy frames are designed heavy enough not 
to deform at the general instability collapse 
pressure. This combination failure of frames 
and shell is not as well understood or formulated 
mathematically as the other two modes, especially 
in the inelastic region. Thus, larger safety 
factors are employed when checking hull designs 
in this mode. 

It should be emphasized at the outset that the program 
to be described designs the stiffened cylindrical shells to 
fail in the axisymmetric yield failure mode. For the 
untested inelastic portion of the program, the design failure 
mode would be axisymmetric inelastic failure. This 
philosophy is used today in submarine design and is 
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•advocated by most naval archietects. It has two basic 
reasons : 

3.) Failu.re by yielding of the shell utilises full 

material (i.e., yield) strength. Stresses in 

buckling type failures are usually below yield 

stress, with failure depending mainly on the value 

of Young's modulus and geometry of the failing 

structure. More efficient structures should thus 

result from designing to a yielding failure, 

2) Imperfections in construction (e.g, , hull out-of- 

roundness) effect buckling failures far more 

2 

seriously than yield failures. Thus, hulls 
designed for yield failure would have less 
stringent requirements for building, and less 
chance of failure, given that imperfections 
might exist. 

The optimization criterion used is a simple 
hull weight to hull displacement ratio, which 
appears to be the best general measure of design 
efficiency for this type of structure. 



Figure 1. 
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FAILURE MODES: STIFFENED CYLINDRICAL 





COLLAPSE FORMULA : REYNOLDS ( DTHB#1 392 ) 




GENERAL INSTABILITY (CONCURRENT FAILURE 
OF SHELL AND FRAMES) 

COLLAPSE FORMULA: KRENZKE~KIERNAN'(DTKB# 16??) 
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TABLE 1. 

LIST OF PROGRAMS USED IN TEE OPTIMIZATION 
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CHAPTER TWO 
HULL STRESS ES 

Numerous subprograms throughout the optimization require 
the values of various stresses in the hull. These stresses 
include circumferential and longitudinal stresses, at the 
frames and at midbay, at the inner and outer shell surfaces. 
Several solutions for stresses due to external hydrostatic 
pressure on ring-stif f ened cylinders have appeared in the 
past. The most famous was due to the Germans, Von Sanden 
and Gunther, in 1920 , Portions of their analysis are still 
in use today. In 1930 » the Italian Viterbo modified their 
analysis to include the so-called stif f ener-expanslon effect 
(a result of axial stresses in the shell). Neither of 
these early analyses, however, included the "beam column" 
effect. This effect, introduced by Salerno and Pulcs^ in 
1951* is caused by the interaction of longitudinal bending 
and longitudinal compression in the hull caused by the axial 
portion of the hydrostatic pressure acting on it. The 
Salerno-Pulos stresses are an exact solution, and the beam- 
column effect accounts for any non-linearities between, 
pressure and strain in the cylinders. In all cases except 
one (see chapter four), the program uses the more accurate 
Salerno-Pulos (hereafter S-P) stresses. 



The beam-column effect is represented in the S-? 
analysis by the parameters (hereafter: GAMA), where: 
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p* is defined'*' ^ as the '’critical load for 
axi symmetric elastic buckling of an unstiffened cylindrical 
shell under the action of uniform axial pressure . 11 GAMA=0 

4 

corresponds to a zero beam-column effect (i,e, , the stress 
solution of Von Sanden and Gunther, hereafter V-G). As 
GAMA grows larger, the beam-column effect, and thus the non- 
linearity between pressure and hull stresses, increases. 

When GAMA > 1, then theoretically the between-f rame failure 
mode shifts to axi symmetric elastic buckling (see chapter 
one for definition). As explained in Chapter 1, this appears 
to be a mathematical phenomenon only, for it has never been 
achieved in actual testing* Other modes of failure (e.g., 
lobar buckling) always appear to occur first, or else the 
axi symmetric failure is always accorcpained by some yielding 
The S-P stresses are calculated by placing various 
combinations of S-P "F functions" (see reference 10 for 
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expressions and curves) into the S~P stress equations. 
Formerly, this was a very laborious process, far more time- 
consuming tha.n obtaining stresses through the simpler V-G 
equations. The curves developed by M.A, Krenzke and R. D. 
Short (included in reference 18) shortened the labor 



considerably. The computerized solution makes it almost 
mandatory to use the superior S-P stress solution. 



The ”F functions” are obtained through the use of 
two basic subprograms (see figures 2, 3 and 4). Subroutine 
FRAME essentially takes shell thickness, unsupported shell 
length, and frame moment of inertia and manipulates them to 
obtain the S-P variables THETA, ALFA, and BETA (see Appendix 
A for method of obtaining frame dimensions). These variable 
are then transferred to .the stress program PULOS via a 
COMMON statement. FRAME is separate from PULOS because in 
one program, (HVYFRM) PULOS is used in two different places 
with the sane scantlings. All variables in FRAME are 
computed in accordance with the S~? stress analysis’*'^, with 
the exception of the expression for effective area, AE. 
Reference 18 lists this as (in program, terms): 

AE-AF*(RM/FF) , or 

AS=AF*(EK/F,F )*-::- 2 , depending on whether 

the framing is internal or external, respectively. This is 

oh 

not strictly correct, and has since been refined by Short' , 
who used the similar equation: 

AE=AF*» ( EM/EF ) , where Q=l+ 2*GMU 

This equation is good for either internally or externally 
framed cylinders with "reasonable” (i.e«, suitable for this 
program's purposes) frame depth/sholl radius ratios. 

The subroutine PULOS is a straightforward adapation 



28 



of the S-P stress analysis. It produces the four "F 
functions” and the variable "a" (see reference 16 ), the 
various combinations of which are used in several parent 
programs to compute the hull stresses. 

From the expression for GAMA, it can be seen that an 
input pressure is required for PULOS to compute its outputs. 
In other words, FRAME and PULOS can be used to calculate 
the hull stresses, given scantlings, material properties, 
and hydrostatic pressure. Lore often, however, the program 
is attempting to f ind a critical failure pressure. This 
results in two unknowns, the pressure and the stresses. 

GALA thus becomes the result of a transcendental process, 
in which PULOS calculates stresses, and the pressure input 
is from a parent program dealing with a particular hull 
failure mode. 

Great difficulties were experienced in this iterative 
process when GAMA^l.O; this meant that ETA1 (see figures 
3 or 4) became imaginary (i.e., the failure mode had shifted 
to axi symmetric elastic). To combat this, several methods 
were employed, each one being different for different 
parent programs. This is the reason for PULOS and PULOS 1. 
The methods of validly circumventing this pitfall will be 
explained separately under each individual failure mode 
subprogram e 
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CHAPTER THREE 



COLLAPSE PRESSURE REDUCTION 

Realistically speaking, no cylinder that is manufacture 
today can be considered "perfect". There will always be 
some reduction in the collapse pressure due to manufacturing 
imperfections, such as shell or frame out-of-roundness , and 
to residual stresses from welding. If, however, a pressure 
hull is inachined rather than rolled and welded, a structure 
that is "perfect" for all practical purposes may be attained 
Of course, for large pressure hulls, the expense (or even 
impossibility) incurred due to size prohibits construction 
by machining only. Thus, come allowances must be made in 
scantling computation. 

The best overall method (i.e,, including perfectly 
plastic and strain hardening materials) yet devised is 
presented as a graph in reference 10 (see figure 5 )* Here, 
the lower curve represents an envelope of numerous model 
test results conducted over the years at the Model Basin. 

It may be observed that "the factor which is all important 
in determining imperfection sensitivity is the margin of 
stability p e /pj| , the abscissa on the plot. The lower the 
margin of stability, the greater the sensitivity to 
imperfections" . The reduction factor (FEDPAC, o1> Pew/pa? 
for machined models ) located on the ordinate is applied to 
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collapse pressures of all modes in the optimization 
program (see figure 6). 

The envelope curve is approximated by the equation: 
REDFAC=0. 64667+0. 11 367*FiT0-0. 0096?*HT0**2 , where 

RT0=p e/ p 1 

This equation and the general form of REDPE were adopted 
from a similar program in reference 22 . 

As Mentioned above, this safety reduction factor is 
applicable not only to this optimization, but also tc one 
for strain hardening metals as well. For this reason (i.e., 
the fact it is easily convertible to the strain hardening 
case) it was chosen above other existing out-of-roundness 
analyses for the strictly elastic buckling cases. No 
mathematical analysis presently exists for the inelastic 

p 1 

case, due to its complexity. - Further discussion on this 
is contained in Chapter 13.. 
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CHAPTER POUR 

A SYMMETRIC (LOBAR) BUCKLING 

The most widely used formulation for this failure mode 

over the years has been the so-called "DTilB Instability 

Formula", developed by Windenburg and Trilling. The formula 

was good only for an elastically-perfectly plastic material. 

Kore recently, Reynolds has developed a more generalized 

formulation which may be used for either elastic-perf ectly 

plastic or strain-hardening materials. In reference 20, 

Reynolds recommended using the V-G stresses (longitudinal 

and circumferential at mid-bay). Ke stated that the accuracy 

of the analysis was not seriously impaired by not using the 

more cumbersome, yet more accurate S-P stresses. This Is 

graphically borne out by the example given in figure 7. 

In any case, the theory correlates very closely with 

A 

experiment, to within four percent . 

The S-P stress programs (FBAKE and PUL03) can easily be 
made common with any other subprogram. Originally it was 
decided that there would be redundancy involved if the less 
accurate V-G stresses were used for the Reynolds lobar 
buckling subprogram (RHLDS). Difficulties, however, were 
immediately encountered when the S-P stresses were utilized. 
RHLDS is used in the main program chiefly as a che cking 
function. Since the hull scantlings, as mentioned in 
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.chapter one, are designed for yield failure, RNLDS is used 
only to insure that the scantling set under scrutiny does 
not fail by lobar buckling. Thus, since scantling sets 
being examined will usually fail first by axisynmetric 
yielding, failure by lobar buckling sometimes occurs at much 
greater pressures. In many cases, failure by axisymmetric 
elastic buckling will occur at a pressure b e tyre on failure by 
yield and failure by lobar buckling. As noted in chapter 
two, this phenomenon causes GAI1A (beam-column effect) to 
increase in value over 1.0, resulting in imaginary values 
occuring within the PUL03 subprograms. Efforts to "force” 
convergence of HELDS by inserting, for instance, values of 
GAKA=0 . 95 , or of taking only absolute values of the radical 
A . 0-GAHA, were only mildly successful. The value of final 
convergence in any case. was not accurate, and certainly was 
not that of lobar collapse pressure. 

For these reasons, then, V-G stresses were used in 
RNLD3 (see figure 8), as originally recommended in reference 
20. One distinct advantage of the V-G stresses' is that 
they require no iteration for convergence. There is no 
separate stress program needed , and the V-G stresses are 
directly (and quickly) computed within HELDS. Since (see 
chapter ten) KaLDS is used itself in an iterative process 
within the main program, this results in substantial savings 
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in computer time. 

Function RilLDS, in this case for elastic-perf ectly 
plastic material, predicts strictly an elastic type faiD.ure. 
Thus, the full pressure reduction from RE DPR (see chapter 
three) is employed directly to compensate for imperfections 
and residual stresses. 



CHAPTER FIVE 



AXI SYMMETRIC (YIELD ) FAILURE 

The real core of the optimization program is the 
axi symmetric failure mode, for the program designs its 
scantlings to fail by yield (see chapter one). Three 
subprograms in addition to FRAME, PULOS, and PiEDPPi are 
included in this grouping: PENCK , ELNCK and THKNS (see 

fig\ires 9f 10, and 11). 

The rather famous Von Sanden and Gunther formulas 92 
and 92 A (utilizing the maximum shell stress theory of 
Rankine at the frame and midbay, respectively) were used in 
design for many years. Recently, however, many more solutions 
have appeared. With the advent of S-? stresses^, the 
stress analysis alone has improved in accuracy. The manner 
in which the stresses are used to predict collapse by 
axi symmetric yield varies greatly. Generally speaking, the 
maximum strain energy theory of Hises and Hencky provides 
the best manner of stress combination to predict failures 
within the various shell yield formulations. The point at 
which this is aoplied is also subject to discussion. 

Although it is generally agreed that the largest stresses 
actually occur at .the frames, it is becoming evident that 
data indicate the best predictors use the mid-bay area as 



the Initiation point of yielding failures. The most 

frequently accurate analysis for ax i symmetric yielding 

1 h 

failure is that due to Lunchik ' ' . Although his formulation 
has not been tested through complete ranges of geometries 
and depths, the tests that have been made indicate his 
solutions are at least as accurate as any others (to within 
l/o of actual failure pressure in many cases), and much better 
than 92 or 92A. in reference 14, Lunchik shows that very 
successful correlations were obtained with tests of ring- 
stiffened cylinders ranging from A (thinness ratio )=0. 4l to 
A =0,70. He recommended his formulation, however, only for 
"cylinders where geometries arc in the range of axi symmetric 
yielding", precisely the case in this program. Basically, 
Lunchik assumes a standard three-hinge failure mechanism, 
postulating that the frame plastic hinges fail first, and 
predicting the pressure at which the mid-bay hinge is 
complete. The basic difference between Lunchik' s analysis 
and others is his computation of "plastic reserve strength". 
His structure does not fail when some outer hull fiber at 
mid-bay has reached yield stress. It fails only after. this 
plasticity has progressed through the shell at that point- 
far enough to produce a hinge and precipitate failure. 

FUNCTION The use of Function RFDP.fi (see chapter 



three) requires both an elastic and an inelastic collapse 
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pressure for the particular failure node under examination. 
For the elastic-perfectly plastic materials, there is no 
"inelastic failure pressure" as such, since no strain- 
hardening is involved. Thus, the yiel d failure pressure 
outlined above is substituted. For the solution of the 
elastic axi symmetric failure pressure, there remain two 
possibilities. 07ie is the "exact" solution offered by the 
S-P stress analysis. When GAMA >1.0 the elastic ax i symmetric 
mode occurs. However, the solution for this pressure is 
bound up in a new "F function" and reqxxires a rather 
involved iteration. Because the REDPR method is approximate 
in any case, an exact elastic axisyrametric value is not 

required. Thus, the S-P solution was rejected in favor of 

1 S 

Lxmchik's inelastic axi symmetric analysis - . This analysis 
is very similar to Reynolds' analysis, for asymmetric 
buckling (see chapter four), in that it can be applied in 
the strain hardening (inelastic case) using the secant and 
tangent moduli. By setting ETAW=0 and ES2C=E, the solution 
breaks down to one for e lastic axi symmetric buckling. This 
is the pressure computed by Function PLHCK, and the process 
is taken directly from reference 13 « 

SUBBOUTIHE KLuC.k The "inelastic" pressure used in 
RSDPfi (and also the pressure to be reduced itself) is 
computed in this subroutine. As explained above, Lunchik's 
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yield analysis for elastic-perf ectly plastic materials 
is used. S-P stresses are used in the formulation (see 
figure 10). This, of course, meant an iterative process, 
and occasionally convergence problems were encountered. 

If input geometries were satisfactory, convergence was 
accomplished in four or five cycles. Occasionally, the 
Iterative program (THKh'S) calling ELNCIl jumped outside the 
range of convergent geometries in its search process. For 
such cases, an iteration limit of 10 was put into ELNCK. 

This meant that the pressure going back to THKNS was not 
entirely accurate , but good enough to continue with a search 
pattern to find a. convergent geometry. An additional 
"safety valve" was built into Subroutine PULOS (see chapter 
two) to prevent GAMA ^1.0 and thus producing imaginary values. 
If GAMA "^1,0, GAMA was set equal to 0,95. and the shell 
thick ness adjusted to achieve this. Thus, depending on 
entering geometries, occasionally shell thickness itself is 
adjusted within ELNCK. This is justifiable, in that the 
overall program is designing to a yield, not a buckling 
failure. Any buckling geometry, even if it could be 
converged upon, would not be desired. Lunch! k 1 s "yield 
pressure", PY (reference 16, equation 35) I s the result of 
the S-P stresses obtained, and Is the pressure at which 
yielding begins at the outer hull fiber, at mid-bay. His 



'plastic reserve strength ratio, PCLEl/PY ( FCTR ) is then 
computed and multiplied by PY to obtain Lunchik ’ s final 
collap.se pressure, PCLE1. The arrival at certain of 
Lunch ik's equations, made somewhat confusing by a misprint 
in reference lL, is done in detail in Appendix B. Before 
sending the collapse pressure to THKNS , it is reduced for 
residual stresses and manufacturing defects by HE DPR . 

At this point (see figure 3 . 2 ), it is interesting to 
see, at least in one case, how the Lunchik and Reynolds 
analyses compare with the log~log plot of hoop stress vs, 
the Windenburg-Trilling formula presented in reference 8. 

For the particular hull diameter chosen, Reynolds’ pressures 
follow Windenburg’s almost exactly. Lunchik ’s pressures 
show hoop stress, at least in this case, to be rather 
conservative. 

SUBROUTINE THRU S This subroutine uses ELNCK in an 
.iterative process to converge on an exact shell thickness 
which will fail by axisyrarnetric yield at the desired 
collapse pressure. If THNLS cannot converge on a thickness, 
it is obvious that the input scantlings, despite changes 
in thickness, are such that failure by elastic axi symmetric 
buckling occurs before axi symmetric yi eld failure. Since 
the full strength of this metal is not then being used, 
this type of solution obviously is 



not desired 



In such 
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a case, THXNS transfers control to the next Iteration loop 
of the main program (see chapter ten). 
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CHAPTER SIX 

HULL THICKNESS IN CREA SE 

Generally speaking, the safety factor introduced by 
REDPB (see chapter three) to 'compensate for "imperfections 
and residual stresses" will be adequate. However, it should 
be remembered that the collapse pressures developed thus 
far all ‘are "triggered" by stress values at midbay. The 
highest stresses actually encountered usually occur at the 
frame faying "flange". To account for this, stresses in that 
area are limited to 755 of yield stress at operat ing depth 
(assumed here to be 2/3 collapse depth) by subroutine 
STRTHK, as recommended in reference 10, This is done to 
account for such things as low-cycle fatigue, creep, stress 
corrosion, and to insure a reasonable stress level in the 

frame flange prior to collapse for those frames with an 

10 

initial out-of -roundness 

Subroutine STRTHK (see figure 13) uses the S-P stress 
1 R 

analysis " to compute all four stresses of interest in the 
shell at the frames: inner and outer plate surfaces, and 

longitudinal and circumferential stress directions. If the 
largest of these is greater than 755 yield stress at 2/3 
collapse pressure, shell thickness is increased in 
increments of 55 until the criterion is met. 
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CHAPTER SEVER 
GE H SRAL 1 i'iST A Bl LITY 

The problem of finding accurate collapse pressures in 
the general instability mode has generally given analysts 
more trouble than the first two failure nodes. The most 
accurate elastic general instability failure analysis was 
done by S. Kendrick in 1953 at Britain's Naval Construction 

i 

Research Establishment, and is generally known as the 

3 

"Kendrick Part III" solu.tion . One year later, A. R. Bryant, 
working in the same establishment , developed a far simpler 

n 

approach/ (Kendrick's was exceedingly complex). Bryant's 

solution was a single two-tern equation which could be (and 

has been) used for design studies without extensive 

computerization, although its accuracy left something to be 

desired (it was non-conservative). Basically speaking, 

Bryant's equation incorporates the "split rigidities 

concept", where one terra represents the contribution of the 

shell, and the other the contribution of the frames and a 

frame-length of shell plating. Although Kendrick's solution 

was put in a simple graphical form by Reynolds^ (later 

3 

extended by Ball"; , it was rejected for use in this program 
for two rea sons: 

1) since the general instability pressure is used 
merely to cjieclc the solution designed for yield 



5 ^ 



failure, extreme accuracy was not required; and 
2) Kendrick’s solution is good only for elastic 
failure, and cannot be applied to strain- 
hardening materials. Since this program is 
designed to be easily converted to one which 
can also handle strain-hardening materials, 
Kendrick Part III was rejected. 

Bryant's solution was also rejected for the second reason 

above. More recently, a very convenient and more accurate 

11 

formula has been developed by Krenzke and Kiernan ". It is 
very similar in structure to Bryant's formula, but can be 
used for either ideally plastic or strain-hardening materials. 
At about the same time, a very similar solution was worked 
out by Lunchlk at the Model Basin. Both the Krenzke-Kiernan 
and the Lunchik analyses give about the sane results when 
compared, with actual test data^. Krenzke-Kiernan was 
somewhat arbitrarily selected for this program, solely 
because it appeared to be referenced more often in the 
literature. 

SUBPOUTl ME KiTKZK This subroutine is called by the 
general instability subprogram, and essentially computes 
failure pressure using the Krenzke-Kiernan formula. It 
also returns the value of "Ij" (number of circumferential 
collapse lobes) which gives the lowest (most critical) 
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failure pressure. It might be noted that this program, like 
some others in the optimization, is designed for dual usage: 
by ideally plastic or by strain-hardening materials (see 
figure lb ). 



SUBROUTINE Gli.'ST (sec figure 15 ) 



tests the 



incoming scantling set from the main program for failure by 
general instability. Depending on whether the collapse 
pressure* is too shallow (test failed) or too d-cep (design, 
too conservative), heavy frame (i.e., bulkhead) spacing is 
adjusted so as to give a collapse pressure that is either: 

1) greater than 1.05 times desired collapse depth 
(this requirement for a small safety margin is 
d\ie to the uncertainties of general instability 
design, and comes from reference 10), or 

2) less than 2 tines desired collapse depth. 

If, in this process, bulkhead spacing becomes longer than 
the entire pressure hull itself, BS is set equal to the 
hull length (BULMTH) , and the requirement for heavy frames 
is dropped (see chapter eight on Subroutine HVYFRh). 
Incoming bulkhead spacing may either be declared or 
undeclared in the input data. If und.e Glared , starting 



bulkhead spacing is taken as approximately twice the hull 
diameter, but always a multiple of the small frame spacing. 
The frame term in the f renzke-Kiernan formula contains, 
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as one of its variables, the "effective moment of inertia", 

El. This is defined as the combined moment of inertia of a 
small frame and its effective lengt h, El,, of shell plating. 
This differs slightly from the Bryant formula, which merely 

uses a frame space (FS) of shell plating. The EL is 

1 P 

determined via the S-P stress analysis- , using the equation: 

EL=SL-::-Fl+B 

The* use of either EL or FS in the formula makes little 
difference in most cases; however, since all the tools were 
handy, the SI, was used when possible. Problems of 
convergence, however, again developed in the PULOS subprogram. 
When general instability pressure became very large (and it 
does for many geometries), the value of GAMA (see chapter 
two) exceeded 1.0. Thus, a separate subprogram (PUL0S1) 
was added, which set EL=FS when GAHA^l.O. This is totally 
acceptable, since this case occu.rs only when the general 
instability collapse pressure is far too deep to worry about 
(its accuracy is only slightly degraded anyway). When 
successive iterations of Gl'NST expand BS and bring this 
collapse pressure- to shallower depths, EL can again be 
accurately computed by PUL0S1 . 



61 



CHAPTER EIGHT 
HEAVY FPAHES 



Subroutine HVYFRH (see figure 1’6) computes the minimum 
size of heavy frames needed to insure the general instability 
collapse pressure calculated by its proceeding subprogram 
GINST. Until recently, heavy frames were designed using the 
standard Levy formula. It wa s found after considerable 



testing, however, that this formula often gave unsafe 

1 7 

estimates . A new formula was derived at the Model Basin 
by Blumenberg in 1965 , which agrees much more closely with 
tests. Results from testing indicated that the effectiveness 
of a particular size of heavy frames decreases as the 
cylinder is lengthened and also that the minimum size of 
heavy frames needed to localize failure between heavy frames 
is possibly not dependent upon their spacing^. Although 
adequate testing has not yet been published to positively 
confirm Blumenberg 1 s formula, initial results show it is 
better than what ias formerly used, and thus it was put into 



HVYFJ i: : 



FIH= 



PCG» BS*RC* *3 

‘( ( * 0 ) -:*ETAI3 j ’ 



where : 



1) E may be substituted for ETAK with ideally plastic 
materials. ETA!: actually is not in Blumenberg’ s 
formula, but was placed there by the author to 
make HVYFFH useful in an inelastic optimization. 
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2) M is the critical buckling node for elastic general 
instability failure of the cylinder with the heavy 
frames replaced by typical frames (i.e., a cylinder 
equal in length, but with no heavy frames). 

(see "List of Symbols" for other variables) 

In computing P.C, the radius to the combined centroid of the 
heavy frame and its effective length of shell plating, an 
ELH had to be determined. This was also provided in 
Blumenberg * s report: 

ELH= SL»F1* (AP-i-SL - *T ) -pry-T 
AFK+SL*T 

It can be seen to be a form of simple ratioing of the 
original EL formula developed in the S~? analysis^. 

Another feature of KVYFRii is to send a value of zero back to 
the main program if the heavy frame spacing computed by 
GINST is equal to (or greater than) the hull length. This, 
of course, would mean that heavy frames are not required, and 
that the hull end closures are used in their piece. 
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CHAPTER NINE 

WSIGHT-DI S PLACE KEN T RA TIO 

The simple weight/displacement ratio was selected for 
the optimization orlterion. Not only was it the simplest 
to use, but it also seemed to be the most general, all- 
encompassing determinant of an optimum design. Due to the 
subprogram system utilized in the optimization, another 
form of criterion could easily be substituted if the need 
arose . 

Essentially (see figure 17)? WTD3P computes the weight 
of a hull section, the length of which is equal to one 
heavy frame spacing, splitting the heavy frame on each end 
in half. This weight is divided by the same hull’s 
displacement in sea water, taking into account internal or 
external frames. 
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DM*T*3S ♦ VF ♦ 
VFH ) 

WTOSP * WT/OISP 



figure 17 
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CHAPTER TEN 
MAIN PROOF AN 



Essentially, the main program performs very few 
calculations. Its primary task is manipulation of the 
various subprograms and iteration control of the entire 
process so as to arrive at an optimum design. See figure 23 
for a basic, easier-to-follou main flow diagram. Figure 24 
is the complete main program flow diagram. In addition, it 
computes "reasonable" entering scantlings for the iteration 
process, based on input data and proven design practice. 

The objective of this type of main program design v;as to 
obtain a system in which the various modes of failure and 
safety factors could be changed or interchanged easily as 
desired. As all failure modes are separated into 
subprograms, it is possible even to substitute an entirely 
different analysis for any individual failure mode. Safety 
factors, 'which so often are subject to modification due to 
new test data, are also easily replaced or changed. 

The input data required (see Appendix D for format) 
are: overall hull length, bulkhead spacing (may or may not 

be given; if not, it computes an optimum B3), internal or 
external framing desired, Young's modulus* metal yield stress, 



metal density, and Poisson’s ratio. 



In addition, as presently 



set up, the program may be iterated for various input values 



6? 



of depth and/or hull diameter. 

The depth input (in feet) is converted to n collapse 
pressure PC (in psi) requirement by using the equation of the 
mean line drawn in figure 18; this represents an average 
depth vs. pal curve for the various oceans included. The 
figure (minus the mean line) was taken from the Handbook of 
Ocean Eng inee ring Tabl es, published by the U. 3. Naval 
Oceanographic Office, and compiled by E. L. Bialek. 

The frame constant CC is computed from the equation of 
the curve shown in figure 19* CC is used to determine frame 
proportions to be used for various depths. The curve is an 
average of proportions of many ring stiffeners used in 
present generation submersible s. The method of obtaining 
frame proportions is printed in the program output (see 
Appendix D). A more detailed explanation may be found in 
Appendix A. 

The program next computes its "datum point", or 
starting set of hull scantlings. Shell thickness, T, is 
computed from the simple hoop stress formula: 

T= PC * 1)11 
*2 . 0*3iGT 

Frame spacing, or, more accurately in this case, 
unsupported length of shell between frames (SL) is confuted 

A A ' L_ 

from the oarameter 0, where: ~ * prrr guvr ■==. 

0= V3( I- T-) 
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Reference 6 lists the usual range of 0 fcr present day 
submersibles as 1.0-2, 5. The main program (see below) is 
set up so that it starts at a naxnium frame spacing and 
then decreases it in further iterations. Thus, a high 
starting value of 0 would be desirable . After several 
calculations at various depths and geometries, a value of G 
=5.0 was used for the datum point. This gives a wide range 
of SL values. The iterative process of the program produces 
values of 0 that eventually go low enough to bracket the 
above range. Other methods of obtaining a starting SL were 
investigated, such as the combined solution of hoop stress 
and the V/indenberg-Trilling equation, but all gave too wide a 
range (usually too large a frame spacing) as depth increased. 
Thus 0 , villi ch does not vary greatly, was selected. 

A starting frame size is selected based on the ratio of- 
frame cross-sectional area to the cross-sectional area of 
one unsupported length of shell plating. Reference 6 lists 
the normal range of this ratio to be 0.2-0. 8. During 
rather extensive investigation, however, it was found that a 
ratio of greater than 0.5 gave frames that were grossly 
overdesigned. Thus, the starting ratio was set at 0.5: 

AF=0.5*SL*T 

The frame moment of inertia is then computed in accordance 
with relations developed in Appendix A. Originally, it 
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x-;a s intended to select a starting frame size by a more 
accurate approach. However » the only formula which could 
be solved in anything approaching a closed form was far too 
conservative (i.e., Tokugawa’s formula, reference 25). The 
advantage of a computerized approach, i.e., investigation 
of a wide range of variables at great speed and low cost, 
was the deciding factor in using the more random iterative 
method described above. 

At this point the program goes into a double loop, 
iterating on SL (outer loop) and AF (inner loop). The 
scantlings are tested by RHLDS. If they do not fail in the 
lobar buckling mode at design collapse depth* the program 
goes on to THaHS. If the scantlings fail BbLDS, SL is 
decreased by multiplying it by 0.9. and the RKLDS test is 
rerun. The decision to run RHLD3 at this point was made in 
order to start as near the ^shoulder" (or at least, not to 
the right of it) of figure 12 as possible, since this is 
generally acknowledged to be an area of optimum design. 
Also, the decision to vary SL, rather than T, to achieve a 
set of scantlings that woxild not fail the RFLD3 test, was 
made for two reasons: 

1) From figure 20, it can be seen that either 
T or SL could be successfully used to change 
J1 j.LD 3 failure pressure, with T being 
slightly more effective. 



RELATIVE EFFECTS OF T*SL OM REYNOLDS' 
LOBAR BUCKLING COLLAPSE PRESSURE AT 
TWO DEPTHS 
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2) The next subprogram in line is THKKS. It is 
connected to RKLDS via another iteration loop. 

To have both programs converging on different 
values of T would cause endless loops in almost 
every case. Thus ; RHJ.D3 was iterated on SL. 

The decision not to vary AP to change RKLDS failure pressure 
is justified by figure 21. The Nott buckling equation is 
almost identical to RKLDS, and thus can be considered 
the same for these qualitative investigations. It may be 
noted that for some geometries, AF has no effect at all on 
PCR; thus, A? was not used in the RKLDS loop. This same 
reasoning (see figure 22) was used in deciding against using 
an AF variation of any kind for convergence within the 
TURKS subprogram. 

As stated a.bove, once the scantling set (revised or 
unrevised) gets successfully past RKLDS, it goes on to 
THKNS. There, the hull thickness, T, is adjusted so as to 
have the shell fail by axl symmetric yield exactly (within 1,1) 
at design collapse pressure. At that point, the scantlings 
with revised T are again looped back through RKLDS to insure 
against lobar buckling, and SL adjusted again a.s necessary. 

If SL has to be re-adjusted, the scantlings again go through 
THiCNS. The process continues until either: 

1) a set of the same scantlings pass without 
change through both RKLDS and THKKS, or 



EFFECT OF FRAME SIZE ON AST 
BUCKLING PRESSURE BY MOTT . 



MME 
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FAILURE PRESSURE (pSi) 



EFFECT OF FRAME SIZE ON AX (SYMMETRIC 75 
(YIELD) FAILURE PRESSURE BY LUNCIIIK. . . 




Figure 22 
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2) SL becomes so snail that there is less than 
four inches clearance between adjacent frame 
flanges. 

In the latter case, the loop is skipped without print-out, 
and the next iteration is started, much the same as what 
occurs if THICKS cannot internally converge on a shell 
thickness (see chapter six). 

i 

Once the scantling set gets past THICKS, T is again 
adjusted upward, if necessary, by STRTHK. A loop to re-test 
back through RMLDS and THKHS is not used here, and would 
serve no useful purpose, since T*s from THKhS and STHTHK 
could rarely ever be made to converge (:l.e., the largest T 
from both programs. is used). 

Prior to entering GINST, the integer 3^ is set equal to 
0. This, used as one of the entering arguments for GIWST, 
is a control variable. 3>0 insures normal (iterative B3 
design) operation of GIM3T. 3>1 is usely solely to obtain 
a single-pass value of H for HVXPRM (see chapter eight). 

The scantling set then enters and is tested by GIKST. 

BS is adjusted as necessary to insure PCG falls between 
1.05 PC and 2PC, the only exception being if BS becomes 
greater than total hull length (see chapter seven). The 
reason that B3, rather than T, SL, or AF, is used here to 
converge on a satisfactory PCG, is because in virtually any 
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'cage, regardless of depth or geometrical proportions, 33 
has a much smaller effect on hull weight for the same PCG 
reduction or increase. Since hull weight already had been 
optimized for shell yield failure (the desired situation), 
any further f, SL, or AF changes would most probably take the 
design seriously off the optimum. 

Once past GlilST, the scantlings pass through HVYFRL 
to obtain a suitable heavy frarae design and WTDSP to 
calculate that particular scantling set’s relative efficiency 
ratio. 

At this point, values of the scantlings are printed out. 
Generally speaking, the print-out of each line of scantlings 
can be said to be optimized on thickness of the shell, 
although certainly other factors (3L, BS , etc.) change as 
necessary to keep the design on a yield failure basis. 

Cnee out of the first loop cycle, AF is decreased by 
•multiplying by 0.8. This is done for ten cycles, so that 
the original AF is reduced to 0.108*AF in the last cycle. 

SL is reduced in the same manner in the outer loop, and a.s 
described above, also may be reduced v.’ithin each iteration 
as necessary to pass BiiLDS. At the end of both loops (100 
major iterations) the program terminates for that set of input 
parameters, and prints out the optimum (i.e. , smallest) 
v,’ eight/displace n e n t r a t i o . 
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CHAPTER ELEVEN 
R ESULTS AMD CONCLUSIO NS 

Various parametric studies were performed Kith the 
program, once it was checked against some contemporary 
submersible hulls to see if it indeed was "in the ball 
park". The most obvious study was to see how the W/D 
ratio varied, with depth (figure 25 ). This was done for 
three common steels in use today, all with approximately 
"ideally plastic" stress-strain characteristics. The points 
obtained plotted into smooth curves on the semi -log plot 
used. It is fairly obvious that if a W/D greater than 0.5 
is considered unsatisfactory, the following limits would 
apply : 



STEEL 


• LIMIT FOR W/D 


HTS 


3500 FT 


HY-80 


6200 FT 


HY-100 


7400 FT 



A second, study was conducted to observe how the W/D 
ratio varied with hull diameter, all other factors 
(including depth) being constant. It was found that it 
varied very little if at all, as can be seen from the 
following set of data; 



WEIGHT/ DISPLACEMENT RATIO 




I 2 3 4 5 10 

DEPTH, THOUSANDS OF FEET 

Fircure 25 



4/To 
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Given: HY 80 steel, collapse depth=5000 ft.: 



HULL DIAMSTER (FT. ) 
20 
30 
50 
80 



W/D RATIO 
0.417 
0,418 
0.417 
0.409 



It could tentatively be said then, that for a given depth, 
hull efficiency (W/D) will be approximately the same 



regardless of how large the hull diameter may be. 

Figure 26 shows the result of plotting the variance of 
W/D with metal yield stress. This also plotted into a 
smooth curve in this case. It is probable, however, that 
for higher yield strength (and therefore, strain-hardening) 
metals requiring a slightly different analysis , the curve 



would have a sharp break. 

The third plot attempted (figure 27) at first appeared 
to be a hopeless scattering of data, but after some analysis 
revealed rather interesting results. The points plotted 
were taken from a single optimization (one diameter, one 
depth) such as in the example print-out in Appendix D. A 
card computing thinness factor ms inserted in the program 



and printed out with the regular data. 

The ”SL loop” is actually a series of W/D points 
computed for scantling sets with the s ame ^ fra me s paci ng , 



WEIGHT/ DISPLACEMENT RATIO 



OPTIMUM V/t/A RATIO VS. STEEL STRENGTH 



05 



CONSTANTS • 

COLLAPSE DEPTH- IOOO FT 
HULL DIAMETER - 30 FT 



Hull length; 2oo FT 
FRAMES •• INTERNAL 




Figure 2 6 



N'GIWLEX 

-4-/TO 



Figure 27 
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buy varying frame size as an input. In each case, it is 
strongly evident that a saddle point, or optimum W/D for 
that frame spacing, was reached. This means that for at 
least this particular set of inputs, the program's frame 
size (AP) iteration range was large enough to bracket 
optimum values. One of these frame spacing minimurns, then, 
wa s the optimum W/D ratio. 

The "W/D envelope" encloses all W/D values computed in 
the program. There seems to be good indication that the 
"optimum W/D" indicated is a true optimum, since the lower 
portion of the envelope rises on either side of it. 

Another indication given by the plot is that the 
program gives a much wider variation in W/D ratios with 



larger frame spacings, for varying frame sizes. As frame 
spacing becomes smaller, the W/D ratios produced become more 
"convergent". 

Perhaps the most obvious conclusion from fig\ire 26 is 
that A alone is certainly not an accurate indicator of 
optimum W/D ratio, although it could be utilized (after 
extensive data gathering) to indicate the general area in 
which to design. 

It is interesting to note that average computation • 
time for each optimum computed was less than one-half 
minute. Also, it was interesting to note that the 
optimum solutions contained as scantlings, generally 
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speaking , smaller frames, larger heavy frame spacing, and 
slightly thicker shells than submersibles of present design. 
Evidently, this results in better W/D ratios. 

One of the major advantages of this program and its 
general method of computation and print-out (see Appendix C) 
is that it gives a great variety of alternatives from which 
to choose, .For instance, perhaps the optimum W/D ratio 
for a certain hu.ll configuration and depth contained 
scantlings which gave a very small frame spacing. 
Analytically, this may give a. superior W/D ratio. But 
practically speaking, the cost of fabricating and installing 
a great number of frames may be out of the question. Thus 
(particularly if the submersible is not critically weight- 
limited), the table of printouts may be "browsed" for 
more attractive scantling sets: ones with accep tible W/D 

ratios, but with inherently lower construction costs. It 
should be repeated, here that each of these printed lines 
are not mere random choices. Each line, prior to print-out, 
has already been through one of the main hurdles of the 
optimisation program. The design of shell thickness and 
frame spacing combination to give failure by axi symmetric 
yield, and thus most efficient use of the material’s 
inhe3?ont strength, is completed prior to printing each line 
of the answer table. Thus, the program nay be used not only 
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as a structural optimization, but also is an indi sponsible 
reference for any economic analysis of a particular design. 
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CHAPTER TWELVE 
RBCOMKBNmTIONS 



This Program. 

1 . Ranges o f Varia tion. One of the most common 

questions after a particular run completion was 
•'1 wonder if this is r eally the optimum, or 
didn't I go far enough in frame (or frame spacing) 
variation? 51 It is believed that plots such as 
figure 27 for each run could probably give a 
definitive answer in most cases, Not only can 
one tell if each set of frames ran through an 
optimum, but -it should also usually be possible 
to tell if the 'W/D envelope" passed through its 
optimum. This may be a rather painful way to 
assure one's self that his run covered all the 
territory that was necessary, but at this writing, 
it appears to be the surest way. As an alternative, 
some rather extensive studies could be conducted 
using various depth and other inputs to test range 
validities. Once determined, the program's range 
of AF and SL variation, governed by the indices K 
and J, respectively, could be controlled 
appropriately. From a rather cursory inspection by 
the author , it appeared that in runs conducted 
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in data-gathering for this paper, the ranges of 
J=K=10 were adequate in all but a few questionable 
cases. 

2. Effects of Fabrication Procedures . It is rather 
obvious that without any doubt, the greatest 
shortcoming in the optimization is in its 
allowances for such relative "unknowns" as residual 
welding stresses, low-cycle fatigue, shell or 
frame out-of -roundness , etc. The gross 
approximations made by REDPR and STRTBK certainly 
fall far short of the accuracy of the various 
failure mode analyses. There is no doubt that more 
basic research must be accomplished in this region 
before a completely dependable optimization 
program could be achieved. As it stands now, the 
safety factors built into the program could fall 
into two categories: 

a) Completely safe design, in which all 
scantlings are overdesigned to the extent 
that design "optimization" is almost useless. 

b) more realistic (i.e., lower) safety factors 
based on scanty experimental data, which 

is not universally applicable, and thus 
might be considered unsafe. 
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It is hoped that the safety devices employed in 
this program adhere to a "middle of the road" 
policy. REDPE, it is believed, results in 
somewhat of an overdesign for the two between- 
frame failure modes, but in underdesign for the 
general instability mode. This underdesign is 
hopefully picked up in STHTHK, which, as mentioned 
in Chapter 6, insures "a reasonable stress level 
in the frame flange prior to collapse for those 
frames with an initial out-of-roundness"^. 

The only way, it appears, to resolve these 
questions, is in extensive testing of models with 
deliberate, measurable defects of a] 1 types. The 
most needed data of this type is in the failure 
mode most effected by defects: general instability. 

p 

It is believed by some authors that the out of 
roundness analyses developed thus far are overly 
pessimisitic when applied to full sized submersibles. 

Fi xtu re Pr ograms . 

The most obvious extension of this program is into 
strain-hardening materials. Because of the 
program’s general, characteristics (i.e., a main 
flow control program manipulating subprograms 
which actually do pressure calculations), it can 
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be easily transformed with .very few changes other 
than subprogram additions. 

The need for such a transformation is obvious, 
when one observes hull steels projected for future 
(and in many cases, present) usage. Any steel 
with a yield strength greater than 100,000 psi 
can be considered to be a strain-hardening material. 
Details of the method of transforming the present 
program into one suitable for use with high-strength 
steels may be found in Appendix E. 



APPENDIX A 



FRAME DIMENSIONS 
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APPENDIX A 
PEAKS DIMENSIONS 

TYPICAL FRAMES : 

In order to arrive at reasonable frame cross-sections, 
averages of a great variety of frame scantlings used, in a 
number of submersibles at various depths were obtained (see 
figure 19, chapter 10). The common denominator for all of 
these frames is the web thickness, B. This method was taken 
from Adame hak in ref erence 1 , who did the same thing for 
surface ship frames. 

The ratios used in figure 19 vie re taken from the 
following, derived from the averages computed (sec figure 

28 ) : 

Flange wldth=FW=CC*B, where CC=y of figure 19- 
Flange depth=FD=FW=CC*B 
Flange thickness=l . 7*B 

In addition to the above, the remaining terms in 
figure 28 are defined as follows: 

PC: Distance of frame centroid from plate neutral 

axis 

DN: Distance of combined centroid from plate 

neutral axis 

The various formulas for frame and combined frame-plate 

gyration wore taken 



moments of inertia and for center of 



T- SECTION 

DIMENSION IDENTIFICATION 
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Figure 28 



M c GiUlXY 

M 70 
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from a general development of these in reference 19 . They 
are listed in easily programmed form, and according to 
reference 19 . differences from values of actual standard 
production frames due to fillets, etc., are always less than 
two per cent. 

The combination of these two methods enables the 
computer to easily arrive at any variety of T-stiffener 
characteristics, each of which is proportioned according to 
present submersible design practice. 

After substituting the expressions for I'D, FW, TF and E 
into the equations for T-fraraes in reference 19» the 
following expressions result for typical (i.e., not heavy) 
frames, and are used in the program: 

Frame are a , AF : 

AF=FVJ#TF+FD*B = 2.7*CC*B**2 
Frame m o ment , of i nertia , FI : 

FI = CC0P*B*#4, where: 

CCCP- ( 2 . 89*CC#*4+11 . 34*CC**3+1 3 . 47*00**2- 
+4 . 36 *CC ) / ( 1 2 , 0*CC+2 0.4) 

Effective fr a me -elate moment of inertia (using 
eff ective lenath of shell nl a tlng, El,, developed 
In ref erence 18 ) , . El : 

KI= 0 . 225*CC*B* «4* ( CC-M . ? ) «-*2* ( ( 1 . 07*CC+0 . $62 ) / 

( COM . ? ) +3 . Of- ( 1 , 63 + ( T~1 . ?,',3 ) / ( ( COM . 7 ) *B ) ) **2/ 

( 1 . Of ( 2 . CO# 3* *2 ) / ( EL-:/f ) ) ) 
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PC: 

FOB* ( ( C C+l . 7 ) / l . 2 2 7 - 1 . 0/1 . 1 7 6 ) +T/2 . 0 



Ml 

D:.=FC-::-2 . 7*CC«B**2/ ( 2 . 7*CCtfIj*tt2+BB*T ) 

HEAVY FRAMES ; 

Heavy frames were averaged in the same manner as 



above, although less data exists. The relationships 
obtained did not depend upon depth, and were determined to 
be as follows: 

FDH=17 . Ov.-BH 
FWH=13.0*BR 



TFH=2 . 0#3H 



CC0P=1 330.0 
Heavy fr a me area, AFH: 

AFH=43. 0*BH**2 

Heavy frame moment of inertia, FIH: 

FIH«1330.0«-B**4 

FCH: 

FCH=1^ . 25#BK+T/2 . 0 
DI-JH : 

DHK=FCH-:fAFfi/(AFH-.‘-2LH-::-T) , where 

ELH~ef f ective length of shell plate next to the 
heavy frame (see chapter 8) 

Fill was not needed in the program. 
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APPENDIX B 

LUNCKIK'S PI A STIC HI MGS ANALYSIS 



(Refer to reference 14 for this discussion) 

Lunch ik ’ s final Pc/Py equation is developed through the 

use of his parameters B and K. There are not solvable except 

by assumptions which Lunch ik makes. Two of his assumptions 

(i.e., and ^ 1 / / o Vl , 0 ' = K/k^p) are easily worked out from 

substitutions in identities from his paper. However, one is 

less clear (i.e., < ^ b ?/c4 yi0 ‘ = ) t particularly because Bp 

is misprinted as in reference 14. 

14 

From Lunchik’s definitions in his analysis of a one 
unit square element: 

(1) M 0 =k?p = circiimf erentia'l edge moment 



(2) Njj=K^ph = circumferential compressive membrane 
force • 

< 3 ) zrhu 

(4 ) cf rn = K 0 p - circumferential membrane stress 

Assuming (i.e., approximating) the circumferential 

bending stress to be elastic: 

- Me v= J2 y . il - _ ^_ M pf tc\ 

^ T' y 2 bh 1 (l)b l 

From ( 2 ) and ( 4 ) : 

' N “ (6, 






Using (5) and (6): 

= 6 Mg h - 

h z N0 hKppk 




6& 0 



^ < 



. D. 
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To further re-arrange 
expressions for B* and. 



(7) 




these assumptions to obtain 
to use in program ELNCK : 

: ( 8 ) K/K, - 





Y<t/n _ h* - 

in?, ' / 1 k. 






/ d'btf \ 



(9) 



(7). (8), and (9) are then put into Lunch ik 1 s equations 
( 12 ), ( 13 )i and ( 15 ) tor G l} 0 a and 0^ , and these in turn 
are substituted in his ( 25 ) and (26), which are vised to 
solve for the plastic reserve stength ratio Pc/Py 



(FCTR in program ELNCK), 
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APPENDIX C 
PROGRAM LISTING 
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APPENDIX I) 



SAMPLE INPUT CARDS 



SAMPLE OUTPUT 



FIRST INPUT DATA CARD 
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Figure 29 



SECOND INPUT DATA. CARD 
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THIRD INPUT BATA CARD 
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E.S.MCGlNLEY,! I 



RING-STIFFENED CYLINDRICAL SHELL OPTIMIZATION 
. . . INPUT PARAMETERS. . . 

MATERIAL DENS I TY=0 .2 85 LBS/CU. IN. YIFLD ST R FS S= 1 OOOOO. PSI 
YOUNGS MODULUS = 0.3E 08 PSI POISSONS RATION. 30 



* FRAME 



DEPTM= 



SHFLL 

THICKNESS 

(In*) 



5c. 0 3 50 
4. 5005 
4. 53 7.4 
4.7713 
4. 0395 
3.8703 
3. 6163 
3. 6163 
3. 6163 
3* 5763 
4. 6769 
4. 3938 
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DIMENSIONS MAY BE OBTAINED BY USING THE FRAME 
CONSTANT BIWEB THICKNESS) AS FOLLOWS: 



.. •LIGHT! TYPICAL ) FRAMES... 
FRAME DEPTH* FLANGE WIDTH* 6.6 X B 
FLANGE THICKNESS* 1.7 X B 

...HEAVY FRAMES... 

FRAME DEPTH* 17 X B 
FLANGE WIDTH* 13 X B 
FLANGF THICKNESS* 2 X B 



5000.0 FT. 


D 1 AMETER 


= 30.0 FT. 


FRAMES: 


INTERNAL 


• 


..DESIGN 


OPTIONS. .. 






TYP ICAL 


TYP ICAL 


HEAVY 


heavy 


WEIGHT/ 


FRAME 


FRAME 


FRAME 


FRAME 


DI SPLACEM 


spacing 


CONSTANT* 


SPACING 


CONSTANT* 


RATIO 


( IN. ) 


( IN. ) 


(FT.) 


( IN. ) 





103.09 


3.438 


200, 00 


0.000 


0. 565 


107.73 


3. 075 


200.00 


0.000 


0.501 


107.40 


2. 750 


200. 00 


0. 000 


0.474 


107.11 


2.460 


200.00 


0. 000 


0.43 5 


106.85 


2. 200 


198.16 


2.981 


0.411 


106.62 


1.968 


73. 77 


2. 664 


0, 3°9 


106.41 


1.760 


47.17 


2.463 


0. 378 


106.22 


1. 5 74 


38. 32 


2.354 


0.374 


106.06 


1 0 403 


20.64 


2. 71 5 


0.393 


105.91 


1.2 60 


20.64 


2.084 


0.378 


86.80 


3.075 


200. 00 


0.000 


0. 531 


86.4 7 


2. 750 


200.00 


0.000 


0.484 



1 2 ? 



4 . 196 ? 


86.18 


2.460 


200.00 


0.000 


0.446 


4.1138 


85.02 


2 . 200 


1 94.23 


3.082 


0.4 33 


3.7526 


85.69 


1.968 


79 . 98 


2 . 677 


0 . 401 


3.7247 


85.48 


1.760 


65.73 


2.651 


0 . 390 


3. 5618 


85 . 29 


1. 5 74 


44.41 


2.399 


0.374 


3.5618 


85 . 13 


1.4 08 


30.22 


2 . 285 


0 . 378 


3. 5228 


84 . 9 8 


1,2 60 


23.14 


2.154 


0.376 


3. 5228 


84 . 85 


1 , 1 27 


1 6 . 07 


2.014 


0 . 385 


4. 2901 


69.73 


2.750 


200.00 


0*000 


0 . 505 


4 . 0 7 3 9 


69.44 


2.460 


200.00 


0.000 


0 . 460 


3 . 8 Q 08 


69.18 


2,200 


200. 00 


0.000 


0 . 422 


3 o 7266 


60.94 


1 .9 68 


190.39 


2.956 


0.401 


3. 5705 


6 8 . 74 


1.760 


70.10 


2.652 


0 . 388 


3.4143 


68 . 55 


1,574 


58.67 


2 , 514 


0 . 365 


3. 4143 


68.38 


1.408 


35.89 


2.354 


0 , 368 


3.4143 


68.24 


1. 2 60 


24 . 51 


2 . 198 


0.372 


3 . 41 43 


69.10 


1.127 


18.83 


2 . 074 


0 . 374 


3.5850 


67.98 


1. 008 


13.17 


1 .998 


0.402 


4 . 00 94 


56.04 


2.460 


200. 00 


0.000 


0 . 482 


3.8607 


55.78 


2.200 


200.00 


0.000 


0.444 


3. 6871 


55.55 


1.968 


199.97 


3.068 


0.416 


3-5315 


55.34 


1 . 760 


80,06 


2 . 660 


0 , 396 


3 . 2611 


55 . 16 


1.5 74 


66.27 


2.596 


0.362 


3.3121 


54 . 99 


1 . 4 C 8 


43.36 


2 . 385 


0 . 362 


3 o 4397 


54.84 


1.2 60 


34.22 


2 * 320 


0 . 369 


3. 4536 


54.71 


1.127 


25 . 10 


2 . 153 


0 . 370 


3.4677 


54.59 


1.008 


16.01 


2.020 


0.383 


3.5216 


54.4 8 


0 , 901 


11. 47 


1. 971 


0.404 


3. 6047 


45.06 


2.200 


94.08 


3.014 


0.474 


3.4426 


44 . 83 


1 . 968 


94.08 


2.845 


0.430 


3. 3388 


44.62 


1.760 


86.65 


2 . 684 


0 . 398 


3. 0832 


44 . 44 


1.574 


68. 13 


2.590 


0 . 362 


3 . 2879 


44.27 


1.408 


64.44 


2 . 537 


0 . 364 


3 . 37 68 


44.12 


1.260 


38.70 


2 . 335 


0 , 369 


3. 3905 


43.9 9 


1.127 


31.37 


2.222 


0.365 


3.4822 


43.87 


1.008 


20 . 40 


2.097 


0.380 


3 . 4963 


43.77 


0.901 


16.75 


1. 971 


0.379 


3.5105 


43.67 


0,8 06 


13 . 12 


1.921 


0 . 388 


3 . 22 94 


36 . 26 


1.968 


67.50 


2.849 


0.449 


2.9829 


36.05 


1. 7 60 


67 . 50 


2 . 708 


0.398 


2. 9176 


35.87 


1.574 


67.50 


2.603 


0.367 


3 . 1114 


35 . 70 


1 .409 
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APPENDIX E 



OUTLINE OF METHOD FOR COM VERTING TO AM 
OPTIMIZATION PROGRAM FOR STRAIN -HARDEN I NG MATERIALS 



Strain-hardening metals differ from elastic perfectly 
plastic metals, in that above the yield point, the stress is 



not a constant value as strain increases, but increases 
(usually at progressively slower rates) as strain increases. 
This means that true plastic flow is never achieved; the 
metal continues to retain a modulus of some value, albeit 
smaller than the original constant value. This results in a 
type of combination elastic-plastic buckling failure in 
stiffened cylindrical shells, termed ine l asti c failure. 

It is particularly important to analyze subersibles 
constructed of strain-hardening materials with strain- 

11 1 3 p A 

hardening analyses' ’ J% ' . Normally, a submarine hull is 
designed to have some plastic yielding (i.e., initially 
beginning in the hull adjacent to the frame flanges) 
somewhere between operating and collapse depth; indeed, there 
must be yielding prior to collapse depth in order for the 
hull to crush there. With ideally plastic materials, once 
the yielding point is reached, generally speaking, buckling 
failure is ruled out (or it would have occurred earlier, due 
to hull geometry). With strain-hardening materials, 
however , once th hull begins to plastically deform. 
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buckling failure is not ruled out. In fact, there is a 
strong possibility that the buckling failure (which will 
occur) will happen at a lower pressure than the yield failure 
calculated by an ideally plastic analysis, due to the 
reduced metal modulus. Thus, to use an ideally plastic 
collapse pressure analysis on a strain hardening material 
might give dangerously overoptimistlc failure predictions, 
particularly if the metal has a very high yield point. 

For convenience in analyses , strain hardening materials' 
stress-strain curves are characterized not only by E, but 
also by ( ETAi'I , tangent modulus) and E s (ESEC, secant 

modulus). See figure 32. 

Because with this type of stress-strain curve, the 

moduli are always dependant upon the stress state when 

above yield stress, it is necessary for all strain-hardening 

collapse pressures to be computed using an iterative process. 

It is generally the approach to the solution of this process, 

20 

and an example pressure analysis by Reynolds that will 
comprise the rest of Appendix E. 

Essentially, the critical collapse pressure is obtained 
when the buckling equation (depending upon buckling node, 
references 3.1, 13 . or 20) is solved simultaneously wi th the 
pre-buckling equation (stress intensity as a function of 
pressure) , where: 
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Figure 32 




Figure 33 
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GRAPHICAL DETERMINATION OF INELASTIC 
BUCKLING PRESSURE 
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Stress intensity^ <S\_ =Jd^~d s -6 % d s (see figure 33) 
Either V-G or S-P stress theories could be used to calculate 
the stress intensity, S. V-G is less accurate, but 3-P 
might give convergence problems, depending upon hull 
geometry and depth at which stress is calculated. 

One difficulty, that of finding a nay of describing a 
strain-hardening stress-strain curve with a minimum of input 
data, is solved in reference 23 . in this method, the entire 
curve may be approximated with extreme accuracy by vising only 
four inputs (see figure 3^) : E, d yt d a and dj, . 3y 

manipulating some of the Ronberg-Osgood equations, it is 
possible to obtain and E s , given any value of stress, via 
the following expressions: 

n= 

lo o io(4/d b ) 

S/E t = l+Ojl-285?n(c5/^) n "' 

E/E s = 1*10 . 428 57(d/d«) 

This will be assumed to be the content of a subprogram 
called Subroutine ROMOS (see figure 35) • 

Both the pro-buckling equilibrium equation and the 
buckling equation will be approximated in the region of 
interest with straight lines. This, and some of the following 
methods of determining the ’ intersection of the two equations, 
were patterned generally after similar methods used in 



reference 22 



METHOD OF OBTAINING ROMBERG -OSGOOD 



INPUT PARAMETERS FOR COMPUTING E s AND E t 




i 
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FLOW DIAGRAM : 
SU-BRCUTIUE ROMOS 



Figure 35 
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Subroutine LINE , This subroutine computes EMM, the 
line slope, and BEE, the line stress intercept of the pre- 
buckling equation, using the S-P analysis (see figure 36 ). 

The values given by LI HE will change only with scantling 
changes. This subprogram would constitute the only deviation 
of the main flow program from the perfectly plastic case. 

3’.t would be placed in the main program directly after point 3 
(i.e., prior to the HMLDS call). It would also be used 
within THL1JS whenever T changes. 

The following discussion will involve the solution of 
the asymmetric inelastic buckling node as developed by 
Reynolds in reference 20. The sane general iterative process 
would be utilized in the solution of the other two modes of 
hull failure. 

Function PCRP . See figure 37 • This small subprogram 
merely computes asymmetric inelastic buckling pressure, 
using as inputs ETAN and ESEC (computed by ROMOS), T, 3L, 
and PCRE1 (computed as PCRE by the elastic portion of PCRP ' s 
calling program , RHLDS). The pressure is computed using 
Reynolds' equations outlined in reference 20 (for a 
simplified presentation, however, see reference 12, which 
gives the axisymraetric mode also). 

S ub routine. .KE PT. See figure 38 . This subprogram obtains 
the intersection of the pre-buckling equilibrium equation 
with the lino determined by the two input pressures (PR31 
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FLOW DIAGRAM: SUBROUTINE LINE 



Figure 36 
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FLOW DIAGRAM : 
FUNCTION PCRP 
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r rijpt 

(PCREI ,T 




FLOW DIAGRAM : 
SUBROUTINE RMPT 



Figure 38 
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and PRS2) and stress intensities (S1G1 and 3IG2) in terras 
of the intersection coordinates, S and P. The values EKK 
and BSE, which determine the pre "buckling equation, are 
piped into RliPT via a CCI : ' r CK statement, along with other 
normal program data inputs, including 3IGA and S1GB, Input 
arguments incl\ide SlGl, SIG2 and PRS1, which along with the 
inelastic buckling pressure at point 2, PRS2 (computed by 
RNPT from input PCRE1 ) , describe end points of the buckling 
.equation line approximation. 

Function BNLD5 . This subprogram is the same as HELDS 
used in the ideally plastic case (i.e., computation of elastic 
lobar buckling pressure, PCREl ) , with the addition of a 
programmed iteration (figure 39 « as used in reference 22) for 
the inelastic portion. The first two stresses used in the 
iteration are SIGY and 0. 9-:iSIGY' . When the difference of 
computed inelastic failure pressure and assumed stress at 
which failure will occur in the next iteration becomes less 
than one half of one per cent of the present predicted 
failure pressure, convergence is assumed. The iteration 
sequence may be followed by using the diagram of the inelastic 
portion of RliLDS (figure 40) viith the plot (figure 39)* 

As noted before, identical procedures would be followed 
in the ax 1 symmetric and general instability cases. The main 



program would be unaltered 



with the 



exception, again as 



noted 



ITERATIVE METHOD USED TO CONVERGE 
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ON INELASTIC COLLAPSE PRESSURE . , . 




P 



COMPUTED INELASTIC FAILURE PRESSURES 

ASSUMED STRESS AT WHICH FAILURE WILL 
OCCUR IN THE NEXT ITERATION 



NOTE: METHOD TAKEN FROM REFERENCE 22 



Figure 39 
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(CONTINUED FROM 

Li/i PitGIj r © • ) 




(Figure kO Continued) 
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before, of insertion 
that an optimization 
much longer (perhaps 
plastic case. 



of subprogram L11\ T E. It is obvious 
using. the inelastic analysis would take 
by five or six times) than the ideally 
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